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A  SEQUENTIAL  PROCEDURE  FOR  SELECTING 
THE  MOST  PROBABLE  MULTINOMIAL  EVENT 
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Francis  Marion  College  and  Clemson  University 

ABSTRACT 

This  paper  deals  with  a  sequential  sampling  procedure  for 
selecting  from  a  given  multinomial  distribution  with  K  .cells, 
the  cell  with  the  largest  probability  of  occurrence.  Observa¬ 
tions  being  taken  sequentially  from  the  given  distribution,  the 
sampling  is  terminated  when  the  largest  count  in  any  cell  is 
equal  to  N  or  when  the  difference  between  the  largest  and  the 
next  largest  cell  counts  is  equal  to  r,  where  r  and  N  are  given 
positive  integers.  The  cell  with  the  largest  count  on  the 
termination  of  sampling  is  selected  as  the  most  probable  event. 
The  probability  of  a  correct  selection  and  the  expected  total 
number  of  observations  are  given. 
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1.  Introduction.  Consider  a  multinomial  distribution  with 


k  >_  2  cells  and  the  associated  probabilities  p^,...,p^  where 

.k 


p^  =  1.  Let  Pj]_j  1  ^[2]  i***±  P £ k ]  ^enote  the  ordered  values 


of  p1,...,pk.  Various  sampling  rules  have  been  considered  in  the 
literature  for  selecting  the  most  probable  event,  that  is,  the 


event  associated  with  p^j*  Beehhofer,  et  al.  (1959)  have  con¬ 


sidered  a  fixed  sample  procedure.  Cacoullos  and  Sobel  (1966) 
and  Alam  (1971)  have  considered  sequential  sampling  procedures. 

The  sampling  procedure  of  Cacoullos  and  Sobel  is  called  inverse 
sampling  in  which  the  observations  being  taken  sequentially  from 
the  given  distribution,  the  sampling  is  terminated  when  the  count 
in  any  cell  reaches  a  given  positive  integer  N.  Clearly,  the 
sample  sice  is  bounded  above  by  KN-K+1.  In  Alam's  procedure 
the  sampling  is  terminated  when  the  difference  between  the  largest 
and  the  next  largest  cell  count  is  equal  to  a  given  positive 
integer  r.  Alam's  procedure  seems  to  be  more  intuitive  than 
the  procedure  of  Cacoullos  and  Sobel.  In  this  paper  we  consider 
a  new  samplir  icedure  which  is  a  combination  of  the  sampling 
rules  of  Alam  a  Cacoullos  and  Sobel.  In  this  procedure  the 
sampling  is  terminated  when  the  largest  cell  count  reaches  N  or 
the  difference  between  the  largest  and  the  next  largest  cell 
count  is  equal  to  r  _<  fJ,  whichever  event  occurs  first.  We  call 
the  new  procedure  R.  As  in  the  case  of  inverse  sampling,  the 
sample  size  for  R  is  less  than  or  equal  to  KN-K+1.  Thus  the 
procedure  R  combines  the  intuitive  approach  of  Alam  with  the\:  \ 

bounded  sample  size  approach  of  Cacoullos  and  Sobel.  V..  T*  ■  *" 


A 


in  each  of  the  above  procedures  the  cell  with  the  largest 


count  at  the  termination  of  sampling  is  selected  as  the  most 
probable  event.  However,  in  the  fixed  sample  procedure  of  Bech- 
hofer  a  randomization  is  necessary  for  selecting  the  most  prob¬ 
able  event  in  case  two  or  more  cells  are  tied  for  the  largest 
count. 

We  consider  the  problem  of  selecting  the  most  probable 

event  from  the  classical  approach  of  ranking  and  selection,  as 

follows:  Let  o  =  (p.,  , . . .  ,o„) .  Given  two  numbers  P*  (4  <  P*  <  1) 
s.  1  K  K 

and  9  >  1,  it  is  required  to  determine  values  of  r  and  N  for 
the  procedure  R,  such  that  the  probability  of  a  correct  selec¬ 
tion  (PCS)  is  at  least  as  large  as  P*  for  all  values  of  p  for 
which 

(p[Kl/p(K-l] 5  -  e*  (1*1) 

A  value  of  p  which  minimizes  the  PCS,  given  (1.1),  is  called 
the  least  favorable  configuration  (LFC) .  It  was  shown  by  Kesten 
and  Morse  (1959)  that  the  LFC  for  the  fixed  sample  procedure 
of  Bechhofer  is  given  by 

Pp  ,  =  (9+K-l)-1,  i=l , . . . ,  K-l 

n  .  a,a...  ,,~1  (1.2) 

PfKJ  '  «  (3+^-1) 

Ce.coullos  and  Sobel  have  shown  in  their  paper  that  the  LFC 
for  the  inverse  sampling  procedure  also  is  given  by  (1.2).  For 
K=2  it  is  clear  that  the  LFC  for  the  procedure  R  is  given  by 
(o..</Om  , )  =  3.  3ut  for  K>2  it  is  not  possible  to  determine 


(3) 


LFC  analytically,  since  the  expression  for  the  PCS  is  mathema¬ 
tically  intractable,  as  is  the  case  with  Alam's  procedure. 
However,  it  has  been  found  through  empirical  evaluation  of  the 
PCS  that  the  LFC  is  given  by  (1.2)  for  K=3.  Therefore,  we  con¬ 
jecture  that  this  result  holds  for  all  K>2. 

In  the  next  section  we  consider  the  probability  of  a  cor¬ 
rect  selection  and  the  expected  sample  size  for  the  procedure  R. 
We  also  compare  the  procedure  R  with  the  three  selection  proce¬ 
dures  mentioned  above. 

Gupta  and  Nagel  (1967)  have  considered  the  problem  of  selec 
ting  the  most  probable  multinomial  event  from  another  approach, 
namely  subset  selection,  in  which  it  is  required  to  select  a 
subset  of  K  events  which  includes  the  most  probable  event. 

2.  Main  results.  The  probability  of. a  correct  selection 
and  the  expected  sample  size  for  the  procedure  R,  given  the 
values  of  r  and  N,  are  derived  as  follows:  the  multinomial  prob 
ability  is  given  by 


.  ,  ,  _K 

g  (x^,  ..  n.  hi=1  , 


(2.1) 


,K 


where  }.  .  x.  =  n.  Let  x'  <  x'< — <x '  denote  the  ordered  values 

— i— i.  i  i  —  z —  —  is. 

of  the  components  of  x  =  (x, , . . . ,x„)  and  let  x  s  A  if  either 

IK  n 

x'  =  N  or  x'  -  xi  i  =  r.  For  x  s  A  let  f  (x, , —  ,>-,)  denote 

K  X  K-i  n  n  1  K 


_  .  .  -iv 

the  total  number  or  ways  or  distributing  n  =  aitferent 


balls  into  K  cells  one  at  a  time,  such  that,  at  any  stage  in 
time  except  at  stage  n,  the  largest  cell  count  is  less  than  N 


(4) 


and  the  difference  between  the  largest  and  next  largest  cell 
count  is  less  than  r.  At  stage  n  there  are  x^  balls  in  the 


ith  cell#  i=l,...,K.  We  have 


PCS  =  irir  ^n  fn (xl' * ’ * 'XK}  P[1]***P[K] 


(2.2) 


where  ?_  denotes  the  summation  over  the  set  of  vectors  x  s  A 

r 

for  which  x„  =  x* .  The  expected  sample  size  S  is  given  bv 

S\  £V 


E(S)  -  in=rK+1  xiAnnfn(xl,**"XK) 


(2.3) 


A  recursive  relation  for  computing  f  (x, ,.-.,x  )  is  given  by 

XI  J. 


fn(xl '  ‘  *  • '  XK}  ~gn  (x! '  *  •  * '  XK}  "im=r^ 


v£B_  (x)  fm {yl ' *  *  *  'yKJ  gn-m  (xl“yl '  *  *  *  ,xK“yKf  1 


where  B^fx)  denotes  the  set  of  vectors  y=  (y^  -  •  •  ,yK)  £  for 
which  i=l,...,K.  The  above  relation  is  satisfactory  for 

computation  only  when  K  and  n  are  small.  We  have  used  it  to 
compute  the  values  of  PCS  and  E(S)  for  K=3.  We  have  applied 
Monte  Carlo  method,  using  1000  trials,  to  compute  those  values 
for  K>3. 


For  K= 2  it  is  fairly  easy  to  derive  exact  expressions  for 

the  values  of  PCS  and  E(S),  as  shown  below.  Consider  a  random 

walk  generated  by  a  sequence  of  i.i.d.  random  variables  (Z  } 

for  which  ?{Z  =1}  =  o  and  P(Z  =-1}  =  a  where  p=p,_,  and  a=o r . , , 
n  -  n  -  -  *  1 2 j  *  *  Uj 

corresponding  to  the  sampling  procedure  R .  Let  S  =  Z.. 

Suppose  that  the  samp line  is  stooped  at  stace  n.  Then  S  =  + r 

or  3  =+(2N-n).  We  observe  that  S  =+r  for  r<n<2JJ-r  and 
n  —  n  —  —  — 

?  =+ (2M-n)  for  2M-r<n«2N-i. 


(5) 


Let  A(n,c;  a,-b)  denote  the  number  of  paths  in  the  random 
walk,  leading  to  =  c  that  do  not  meet  the  line  y  =  a  or 

y  =  -b  where  -b  <  c  <  a  is  given  by  (see  e.g.  Billingsley  (1963)  , 
Formula  (11.5),  p.  78) 


m  w 

A(n,c;  a,-b)  =  Ln,c"iK=1(Ln,2  (K-l)  (a-i-b) +2a-c“Ln,  2X  (a+b) +c 


+  Ln,2(K-l) (a+b)+2b+c“Ln,2K(a+b)-c) 


(2.4) 


where  L  =  (  )  denotes  the  total  number  of  paths  leading  to 

n  t  c  n  **c 

2 

=  c.  From  (2.4)  we  have  after  simplification 

n-1  n-1 

A (n-1 , r- 1 ; r , -r )  =  (^J  -  >K=1  ^(4K-3)r4-n) 

2  2 

+  (  n“1  )  -  (  n“1  )  -  (  n“1  )} 
(4K-1)  r-i-n-2  (4K~l)rtn  (4K+1)  r+n-2 

2  2  2 


(2.5) 


,N+r-l1 


A<n-l,29-n-l*r,-r>»  (£*>  -  ^  *r  5  <  (2K-1 )  r-W 


*  (  (2K-l)r^N-l}  "  (2Kr+M-l;  “  (2Xr-N+n}  ;  U’6) 


where  [yj  denotes  the  greatest  integer  loss  than  or  equal  to  y 
and  (^5  =  0  if  y  >  n-  Associating  the  random  walk  with  the  pro¬ 
cedure  R  we  derive  the  following  result: 


(6) 


PCS  =  v 


n+r  n-r 
0  2 

A(n-i,r-l;  r,-r)p  q 


n=r,r+2,...,2N-r  "WJ  x 

2N— 1 

'n_2h  -+1  A(n-l,2N-n-l;r,-r)pN  q 


M  n-N 


E{S)  =  E 


n=r, r+2 , . . . , 2N-r 


1 , 2N-n-l; r ,-r) p  q  (2.7) 

n+r  n-r  n-r  n+r 
nA (n-l,r-l;r ,-r)  (p  2  q  2  +D  2  q  2  ) 


+  Zn=2N-r+l  nA (R-1 ' 2-;‘n_1  ? r » _r)  (p V‘\dN  V'3)  (2.3) 


The  formulas  for  PCS  and  E(S)  are  simplified  for  r=2,3 


and  N-l ,  as  given  below. 

K=2 , r=2. 

™  .  P2(1-(2Pct)N-;1)  .  .N-l  N  N-l 
PCS - p—  +  2  P  q 

n+2  n-2  n-2  n+2 

E(5)  =  I  .  .  .n.2  ,2  "T”  2  2  . 

n—  2  /  4  t .  . .  1 2N— 2  {p  q  +p  cj  ) 


K=2 , r=3. 


xnv  .  .  .  N  2M-1  2N-1  N 

+  ( 2N- 1 )  (p  q  +p  q  )  . 


P3(1-(3p<t|N“2)  .  ,!i-2  MS-2 _ 

- TT55S -  +  3  ?  9  (lT2qJ 


E  (Si  = 


n=3,5,...,2N 


n-3  n+3  n-3  n-3  n+3 
2  “>  2  o  o 

.  n.  3  (p  q  *P  q  “  } 


-  ,  I*  -«  4.  ,  .*  t  W.  _  -,S<  —  —  ,  ii  iH—  i.  It—  J.  0. . 

3  (p  q  +p  q  ) +2 . j  (p  q  +p  q  ) 


.  O  y _ V— ,  1 

A—  »  L—-\—  i  • 


N-l  >j  N-2  N-ri-1  . 

PCS  =  p*  S(N-15?  cr+Z  -_-,((  ,  )-1)?  <j] 

■*  }-4 


+  {  {  zs*'~k'~x 

1  '  N-l'  ~ 


,  N_ i 

(N,N)  -o  *s  "(2T5-1J 


(7) 


M_  M_  I  2  M  M 

E(S)  =  (N-l)  (p  ’+<T  }  +  (N  -1)  (p  q+pq  ) 

+  Z  (N+j)  (N^~1)-l)pN'C!^ 
j=2  M“A 


,  lw/2M-2,  ,,,  N.  N-l 

+  (2N-1)  (  {  ,  >  -2)  p  q 


where 

x  \  r(m+n)  f ®  m-l.n  .n-l, 

Ip(rn,n)  P  (m)  (n)  0  x  (1-X)  dx 

denotes  the  incomplete  beta  function. 

Table  I  gives  the  values  of  PCS  and  E(S)  for  the  least 
favorable  configuration  given  by  (1.2)  and  specified  values 
of  r,N  and  8.  The  values  of  PCS  and  E(S)  are  given  in  that 
order  within  parentheses  in  the  table.  The^  figures  for  K  >  3 
were  computed  from  the  results  of  a  simulation  study  based  on 
1000  repetitions  of  e  sampling  process. 

In  the  classica  formulation  of  the  selection  problem  it 
is  required  to  determine  the  values  of  (r,N)  for  which  PCS  >  ?*, 
aiven  d.v1/o...  >  e,  where  ?*  and  8  are  given  numbers.  For 

*  [Kj  *  f/.-ij  — 

each  r  Twe  select  the  smallest  value  of  N  for  which  PCS  _>  ?*. 

There  may  be  several  pairs  of  values  of  (r,N)  for  which  PCS  >  ?*. 
However,  it  has  been  found  empirically  that  the  pair  corresponding 
to  the  smallest  value  of  r  leads  to  the  smallest  value  of  E(S), 
Therefore,  we  select  this  pair.  The  following  table  gives  values 


of  such  pairs 

(r,N)  and  ECS)  for 

d*=  7 ^  on  £  —  i  £ 

2.0,  2.4, 

j.O  rv=2{i)o. 

The  figures  given 

in  the  table  for  K  > 

4  are  based 

on  Monte  Carle 

>  results. 

Monte  Carlo 


(8) 


Values  of  (r  ,N,E  (S) )  for  PCS  >.  P* 

P*=. 75 


K: 

2 

3 

4 

5 

6 

8=1.6 

(3,5,5.96) 

(4,12,17.19) 

(4,19,30.31) 

(4,23,41.08) 

(5,25,66.34) 

2.0 

(2,3,3. 09) 

(3,6,3.44) 

(3,9,13.56) 

(3,11,13.76) 

(3,16,23.84) 

2.4 

(2,2,2.42) 

(2,  5,4.70) 

(3,5,9.07) 

(3,6,12.57) 

(3,6,14.53) 

3.0 

(1,1,  1) 

(2,3,3.63) 

(2,5,5.58) 

(2,5,6.33) 

(2,3,8.30) 

?*=. 90 

- 

3=1.6 

(5,21, 17.00) 

(6,  30,34.77) 

(10,28,61.95) 

(3,31,84.09) . 

(8,36,112.20) 

2.0 

(4,3,3.90) 

(4,15, 15.01) 

(4,19,23.39) 

(4,21,30.49) 

(5,20,47.87) 

2.4 

(3,6,5.72) 

(3,  20,3.90) 

(4,9,15.91) 

(4,11,21.58) 

(4,13,27.61) 

3.0 

(2,14,3.20) 

(3,5,6.  44) 

(3,8,10.02) ‘ 

(3,3,12.84) 

(3,9,14.80) 

To  compare  the  procedure  R  with  the  three  procedures  men¬ 
tioned  in  the  previous  section,  the  following  table  shows  the 
minimum  expected  sample  size  rounded  to  the  next  largest  integer  fo 
the  four  procedures  corresponding  to  the  least  favorable  config¬ 
uration  given  by  (1.1)  for  9  =  1.6,  2.0,  3.0  for  which  PCS 
P*  =  .75,. 90.  We  denote  in  the  table  the  sample  size  for  the 
fixed  sample  procedure  of  Bechhofer  et.  al.  [2]  by  F  and  the 
expected  sample  size  for  the  sequential  sampling  procedures  of 
Cacoullos  and  Sobel  14]  and  Alam  [1]  by  E { I >  and  E (A) ,  respec¬ 
tively.  The  values  of  F  and  Ed)  are  reproduced  from  Table  II 
of  Alam  [1].  Note  that  the  procedure  R  is  the  same  as  Alam's  pro¬ 
cedure  when  it  is  seen  from  the  table  that 


(10) 
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(11) 


Table  I-Probability  of  correct  selection  and  expected 

sample  size 


K=2 


e  =  1.2 


N 

10 

15 

20 

25 

r-2 

(.590,  4.0) 

(.590,  4.0) 

(.590,  4.0) 

(.590,  4.0) 

4 

(.  647,11.8) 

(.662,13.8) 

(.669,14.7) 

(.672,15.1) 

6 

(.656,15.4) 

(.686,21.0) 

(.706,24.8) 

(.720,27.4) 

8 

(.656,16.3) 

(.689,24.1) 

(.715,30.5) 

(.735,35.6) 

6  -  1.6 

2 

(.719,  3.8) 

(.719,  3.8) 

(.719.  3.8) 

(.719,  3.8) 

4 

(.832,10. 8) 

(.'855,12.15 

(.863,12.5) 

(.866,12.7) 

6 

(.848,14.3) 

(.892,18.4) 

(.916,20.5) 

(.929,21.7) 

8 

(.849,15.4) 

(.898,21.8) 

•(.927,26.0) 

(.946,28.6) 

10 

(.849,15.6) 

(.898,23.3) 

(.929,29.4) 

(.950,33.6) 

■ 

o 

• 

eg 

II 

CD 

N 

6 

9 

12 

15 

<N 

II 

V, 

(.800,  3.6) 

(.  800,  3.  6) 

(.800,  3.6) 

(.800,  3.6) 

3 

(.860,  6.1) 

(.880,  6.7) 

(.886,  6.9) 

(.888,  7.0) 

4 

(.875,  7.7) 

(.912,  9.3) 

(.929,10.0) 

(.936,10.3) 

5 

(.878,  8.4) 

(.922,11.1) 

( .  945,12.5) 

(.957,13.3) 

6 

(.878,  8.6) 

(.924,12.3) 

(.950.14.5) 

(.965,15.8) 

Best  Available  Copy 


(12) 


e  =  2.4 


6  ■ 

9 

12 

15 

gsg- 

1  * 

(.850, 

3.4) 

(.852,  3.4) 

(.852,  3,4) 

(.852,  3.4) 

(.911, 

5.7) 

(.927,  6.2) 

(.931,  6.3) 

(.932,  6.3) 

(.926, 

7.2) 

(.955,  8.5) 

(.965,  8.9) 

(.969,  9.1) 

1 

(.928, 

8.0) 

(.962,10.2) 

(.978,11.2) 

(.983,11.5) 

Sfi:" 

(.928, 

8.3) 

(.964,11.5) 

(.981,-13.0) 

(.988,13.8) 

■ 

w. 

(.899, 

3.2) 

0  =  3.0 

(.899,  3.2) 

(.899,  3.2) 

(.900,  3.2) 

I 

(.952, 

5.3) 

(.962,  5.5) 

(.964,  5.6) 

(.964,  5.6) 

I 

I 

(.964, 

6.7) 

(.982,  7.5) 

(.986.  7.7) 

(.987  ,  7.8) 

■ 

■ 

(.966, 

7.6) 

(.986,  9.2) 

(.993,  9.7) 

(.995,  9.8) 

a 

i 

|  (.966, 

7.9) 

(.987,10,5) 

(.995,11.5) 

(.997,11.8) 

i 

■ 

i 

J 


(15) 


Table  I  (Continued) 


0  =  1.2 


(.360,8.14)  (.360,8.20)  (.360,8.20)  (.360,8.20) 

(.417,23.80)  (.460,31.61)  (.479,35.65)  (.483,37.83) 

(.426,27.39)  (.460,41.85)  (.490,53.87)'  (.509,63.35) 

(.426,27.69)  (.461,44.11)  (.488,59.98)  (.517,74.42) 

(.426,27.70)  (.461,44.44)  (.488,61.28)  (.516,77.85) 


0  =  1.6 

(.523,7.32)  (.521,7.86)  (.523,7.75)  (.501,7.58) 

(.632,21.99)  (.731,27.25)  (.762,30.10)  (.771,30.48) 

(.670,25.16)  (.769,37.06)  (.814,45.28)  (.840,50.80) 

(.679,26.19)  (.773,40.24)  (.847,51.59)  (.878,62.41) 

(.675,25.30)  (.779,41.21)  (.840,54.73)  (.865,66.13) 


3  =  2.0 

(.643,7.17)  (.625,7.28)  (.626,7.01)  (.666,7.35) 

(.765,14.41)  (.806,15.33)  (.798,15.87)  (.798,16.05) 

(.831,19.33)  (.888,22.65)  (.900,23.42)  (.894,23.84) 

(.838,22.05)  (.929,27.66)  (.924,31.14)  (.942,31.00) 

(.854,23.02)  (.924,31.83)  (.951,35.16)  (.965,38.18) 


WM 


(16) 


USsfr-  % 

■g=n; 


8  =  2.4 


10 

15 

20 

25 

(.700,6.48) 

(.723,6.57) 

(.723,6.71) 

(.718,6.48)  • 

(.  865,12.35) 

(.879,13.00) 

(.878,13.73) 

(.871,13.36) 

(.925,17.34) 

(.949,18.77) 

(.945,19.51) 

(.953,18.73) 

(.931,19.77) 

(.966,23.87) 

(.977,23.79) 

(.982,24.99) 

(.930,20.95) 

(.968,26.56) 

(.983,28.43) 

(.985,29.34) 

8  =  3.0 

(.805,5. 81) 

(.808,6.04) 

(.816,5.66) 

(.820,5.79) 

(.917,10.52) 

(.928,10.78) 

(.937,10.86) 

(.935,11.22) 

(.959,13.71) 

(.980,14.80) 

(.969,15.19) 

(.976.14.80) 

(.974,16.55) 

(.992,18.48) 

(.992,19.27) 

(.994,19.08)  \ 

(.966,18.43) 

(.994,21.64) 

(.996,21.91) 

(1.000,22.49)  i 

1 

— — -  I  m 

p_|¥] 


(.290,11.04) 

(.313,10.37) 

(.319,10.88) 

(.309,10.32) 

(.351,29.25) 

(.366,40.54) 

(.371,46.82) 

(.387,52  23) 

(.345,32.50) 

(.366,51.26) 

(.436,66.49) 

(.453,80.13) 

(.357,32.76) 

(.370,53.60) 

(.435,72.40) 

(.410,91.60) 

(.355,33.00) 

(.371,53.20) 

(.432,73.88) 

(.455,95.02) 

8  =  1.6 

(.459,9.79) 

(.439,9.95) 

(.446,10.31) 

(.466,10.38) 

(.599,27.60) 

(.661,34.60) 

(.679,38.86) 

(.734,41.61) 

(.626,30.85) 

(.727,46.06) 

(.798,56.41) 

(.837,63.70) 

(.602,30.96) 

(.749,48.83) 

(.798,63.84) 

(.843,76.92) 

(.646,30.64) 

(.671,49.53) 

(.791,66.64) 

(.844,82.71) 

6  =  2.0 

(.573,9.10) 

(.590,8.93) 

(.562,8.88) 

(.556,9.26) 

(.735,17.68) 

(.777,19.48) 

(.780,20.87) 

(.771,20.14) 

(.775,24.16) 

(.864,28.94) 

(.885,30.15) 

(.870,30.21) 

(.783,26.64) 

(.395,34.90) 

( .  9  30 , 38 . 65)_ 

(.950,40.77) 

(.799,23.12) 

(.904,33.99) 

(.943,45.19) 

(.967,47*41) 

{.636,3.10) 
{.839,15.59) 
{.832,21. 11) 
(.899,23.78) 


9  =  2.4 
(.658,8.04) 

( . 854 ,16.66) 
(.922,23.79) 
(.960,23.75) 


(.672,8.30) 
(.862,16.77) 
( .947,24.19) 
(.954,30.69) 
(.980,35.80) 


(.671,8.25) 

(.864,17.20) 

(.944,25.00) 

(.979,31.09) 

(.985,36.65) 


(.887,25.51)  (.965,32.53) 
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